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■ Abstract 

'■ We present a variational method which uses a quartic exponential function as a trial 

^ I wave-function to describe time- dependent quantum mechanical systems. We introduce a 
if^ I new physical variable y which is appropriate to describe the shape of wave-packet, and 
^ I calculate the effective action as a function of both the dispersion ^ (g^) and y. The effective 
potential successfully describes the transition of the system from the false vacuum to the 
^ ■ true vacuum. The present method well describes the time evolution of the wave-function 
of the system for short period for the quantum roll problem and describes the long-time 
evolution up to 75% accuracy. These are shown in comparison with the direct numerical 
Q^l computations of wave-function. We briefly discuss the large N behavior of the present 
^ ■ approximation. 



PACS number: ll.80.Fv, 73.43.Nq, ll.30.Qc, ll.15.Tk 



Phase transition is one of the most important physical phenomena in nature and has 
wide range of applications to condensed matter physics, particle physics, and cosmology. 
Most of the studies on this subject have been done in the framework of quasi-static 
transition or using the Gaussian ansatz developed by Jackiw and Kerman pQ . There have 
been many attempts 013111 to go beyond the Gaussian approximations. It is our purpose 
in this paper to go beyond the Gaussian approximation in two respects. First, we need 
a fully non-perturbative way which links the initial Gaussian packet (GP, false vacuum) 
to the symmetry broken degenerate vacuum state (true vacuum). Second, we try to find 
the relevant physical parameters which describes the symmetry breaking effectively. 

In this paper, we consider a quantum mechanical model for time-dependent dynamics 
described by the potential, 

vm = ^[f-e{t)]\ (1) 
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where k'^{t) increases from a negative value to a positive number k'^ asymptotically. The 
initial GP centered at g = cannot remain as Gaussian during the time-evolution, but 
evolves to the packet centered around two minima of the potential as fc^(t) approaches 
K^. For ^ oo, the new ground states are linear sum or difference of two un-correlated 
GPs centered at each minima. In this case, the two ground states are degenerated. 

The dispersion (g^) of a wavepacket may describe the size of a GP or the distance 
between two packets of a double Gaussian packet (DGP). To discern the shapes (for 
example, GP or DGP) of wavepackets of the same dispersion we introduce a dimensionless 
quantity y, which we call "shape factor", in addition to the dispersion (g^): 
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(2) 



A similar expectation value as y was calculated in Ref. [3] in relation to the new infla- 
tionary scenario. To illustrate the role of y variable, consider a wavefunction which is a 
sum of two GPs of the same size. If ?/ = 1, the density of each GP is a delta function 
or the two GPs are infinitely far away so that no correlation exist between them, which 
provides the lower bound of y{> 1). If the two GPs completely overlap, it corresponds to 
?/ = 3, a single Gaussian packet. In between these two states, 1 < ?/ < 3, the two GPs are 
mixed and interfere with each other. For ?/ > 3, there are no separable packets, and the 
wavefunctions are better localized than GP [H]. 

The effective action in the variational method PP is given by 



dt{ip,t\idt - H\ip,t), 



(3) 



where H = ^ — \-V{q,t) and we use ^ = 1. In this paper we use the trial wavefunction. 



{Q\ij,t) = N-'exp 



2 V 2/i^ 



(4) 



which has both of the DGP (x oo) and the GP (x — oo) limits, where we assume 
/i > 0. In the static case, the double Gaussian approximation was used in Ref. [7j, where 
a sum of two Gaussian functions is used as a trial wave-function. However, it is difficult 
to generalize the double-Gaussian method to the case for time- dependent systems. The 
normalization factor can be determined by the following integral: 



= 

where f{x) is given by jH| 



dQexp I -^r^ + 



2/i^ 



2 TT 
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I^iix'^) + sgn{x)Ii{x^) 



The dispersion and the "shape factor" for this wavefunction are 



1 + 2x/7/ 
/'V(2/2) ■ 



(5) 



(6) 



(7) 
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y{x) is a non-increasing function of x from 3 to 1, which makes the inverse function x{y) 
be defined uniquely. We use ?/ as a basic variable instead of x, because its range is bounded 
below hj y = 1 for any kinds of wavepacket j6jj and it has definite physical meaning. The 
expectation values of other polynomials of can be written in terms of these parameters. 
With this trial wavefunction the effective action is given by 



dt 



yq'^tl 2 

— q p-2y 
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where Y{y) = 2xf'/f = yf'"^ /{2p) — 1 and the effective potential is 



imq^ 



with the free potential, Vp, given by 

Vpiy) = 1 + 



(3-i/)(y + i) 



(9) 
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This free potential, coming from the expectation value (p^), represents the effect of quan- 
tum mechanical uncertainty. The expectation value of symmetric potential, V{q, t) = 
Voit) + \k{t)q^ + + ^q^ + ■ ■ -, with respect to [V', t) is 



{V) = Vo{t) + 4^g2 + -^yq" + -^y 



4! 
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From the action (jH]), we notice that 11 and p are the momentum conjugates to —yq'^/2 
and g^, respectively. 

Let us solve 11 and p equations first: 



\n{yq^) 



y 



Inq 



y(l + Y) 
yq^H p 
m 



yq'^U 



m 



p_ 
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Removing 11 and p by Eq. (fT^ is just the Legendre transformation. Introducing new 
drj 



variable t] by 
and g. 



dy 



D 



I 1 + Y 
z/(3 - y) 



we get a quite simple effective action in terms of r] 



S 



dt 



mq 



m 



2 ^ 



U[q,y{v)] 
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The dynamical equations of motion for q and rj are given by 



m 



d_ r 
dt 



q^7] 
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D 



mq 



Qmq^ 

Vf 
Amq^ 
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(14) 
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The free potential, has an absolute minimum at = 3, q = oo) and is positive 

definite. An interesting point here is that y = 3, the GP, actually corresponds to x = — oo. 
On the other hand, in the effective potential (^UJ, ?/ = 3 is a regular point, which can 
be extended to larger values. This property of the effective potential implies that the 
trial wavefunction is insufficient to give a full description of y-dependence and we 
need a more general trial wavefunction for complete quantum mechanical description 
which includes the range y > 3. The generalization of the trial function to = 
iV~^(l + \z\Q'^)~^\{Q\'iIj ,t)\'^ may allow a full description of y-dependence, but, is difficult 
to integrate. Instead, we use a patch for the region ?/ > 3, 

{Q\i,,t) = N-^ew^-Q^ + iTl^Q^ + {-x + ip)\Q\Y y > 3, (15) 

where we assume x > and G > 0. The wavefunction is singular at the origin Q = due 
to the IQI term in the exponent. Hence we regard the absolute value as a small a limit 
of y/Q'^ + a?- With this trial wavefunction, we have the same effective action as Eq. (|T!^ 
d tt 

with D = - — ^ , {y > 3), and 

2a\/a — 1 



2 

l + 2z--.. 

/ V TT 




for y>3, (16) 



where z = 2x^G, f{z) = e^(l — erfy^), and a = ^^{^^'^ ^-^. The "shape 

2z[ll[^Tizf) - IJ^ 



factor" y{z) monotonically increases from 3 to 6 as a function of z, 

, , _ 3 + 12^ + - 2(5 + 2z)v/i/(v/?/) 

[1 + 2. - 2 Vi/(v^/)F ■ ^ 

As an example, let us consider the time-evolution of an initial wavepacket given by 
Eq. (jH) with {y = yo,q = qo) in the harmonic potential s{t)q^ /2. The dynamics of y can 
be evaluated exactly from the elliptic integral 



JriO 



where = 4m^[?7(0)g^(0)]^ + Vf[?/(0)] is an integration constant. The equation of motion 
for q becomes 

m = -smt) + (19) 

If the system is potential free [s{t) = 0], ri{t) asymptotically approaches to a fixed value 
rjf and q{t) asymptotically increase with a constant velocity determined by the energy 
conservation law. The allowed range, for constant s{t) = s, of g and y is also determined 
by the energy conservation law: 

^^ + \sq'<E,,, = Uiqo,yo). (20) 
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Let us now consider the effective potential for the classical potential (0). The effective 
potential Q naturally determines the true ground state with the condition, 



d,U{q,y)=0 = dyU{q,y). 



(21) 



To see the behavior of the effective potential more clearly, we variationally determine y, 
and then write down the effective potential in q: 



Veffiq) = U{q,y^{q)) 



VriyM) , A r 
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q'y^{q)-2e{t)q' + k\t) 



(22) 



where yy means that we determine y by minimizing the effective potential with the condi- 

In the static case this value corresponds to the minimum position 



Smq'^ 



tion, ^q^ - 

of the potential for a given q, which takes |/„(0) = 3 and yv{oc) = 1. We present figures 
(Fig. 1) of the effective potential for two typical sets of parameters. 



Fig. la 



Fig. lb 
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Figure 1: The effective potential as a function of q for the parameters A = 0.1, k = 
7, m = 1 (a), and A = 123, k = 1, m = 1 (b). 

In these figures, l^// = U, Vcq, and Vc^g^ represent the effective potential (j22|l . 
the Gaussian approximated effective potential for (g), and the Gaussian approximated 
effective potential for (g^), respectively. Here, V^q can be calculated from Eqs. (2.9) 
and (4.6) of Ref. jS] with slight notational change (g = (Q)) and V^qZ from Eq. (2.9) of 
Ref. |3] with {Q) = and G q^. The effective potential, K//, is very close to Vc^q^ for 
q < k/3, and it becomes close to Vc^q for q > 2k/3. This clearly shows that the initial 
GP is divided into a DGP, with each packet of the DGP moving as if it is a free GP for 

large k. The value of 1 + (]^^|^) , for ~ 1.] is effectively 1 for the most range 
of g if is sufficiently large (Am)~^/^], since the characteristic size of g is 0{k). 
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Let us explicitly describe the dynamics of an initial GP with q = qo for the time- 

dependent potential (1). Because of the transition (as k{t) increase) y eventually goes 
to 1 for most part of the dynamics. The potential energy difference AV^ = Veff{qo,y = 
1) — Veff{qo,y = 3) > and the presence of kinetic energy in y prevent q from reaching 
go- The time-dependence of k{t) decreases the total energy so that q oscillates near the 
true vacuum. We present, in Fig. 2, a solution of the differential equation (fT^ and 
its exact numerical solution for the case of k{t) linearly increasing to a finite value for 
about a half period of q. In this example, we do not need the patching process by the 
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Figure 2: Solution (Left) of Eq. ()14|) and exact numerical time evolution by wave-function 
simulation (Right) of q{t) and y{t). In this figure, we set A = 0.0123, m = 1, go = 1-6, 
and y{0) = 3. k{t) = t/3 during < t < 15 and remains constant afterward. 

wavefunction (jl5|) since states with y > 3 do not appear. To see the long time behavior 
of the system, we present one more figure(Fig. 3). The main characteristic feature of the 
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Figure 3: Solution (Left) of Eq. ()14|) and exact numerical time evolution by wave-function 
simulation (Right) of q{t) and y(t). In this figure, we set A = 0.0123, m = 1, go = 1.6, 
and y{0) = 2.906. k{t) = t/3 during < t < 30 and remains constant afterward. 

long time behavior is the oscillation of the amplitude of short time oscillation. The error 
of the oscillation period in Fig. 3 is about 25%. This error comes from the non-exactness 
of the variational wave-function to the exact time-evolution of the wave-function. 
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As another example, let us consider the quenched transition potential with k{t > 0) = 
K and k{t < 0) = 0. The discussions above for the time-dependent k{t) transition also 
applies to the present example. We present a numerical solution of the differential equa- 
tion (|T^ and its exact time evolution in Fig. 4. The state with large y{> 3) periodically 





Figure 4: Solution (Left) of Eq. (fT^ and exact numerical time evolution (Right) of q(t) 
and y{t) by wave-function simulation. In this figure, we set A = 0.0123, k = 5, m = 1, 
go = 1-6, and y{0) = 3. 

appears. This means that we need the patching process (fTK|) for the time evolution of 
this system. Comparing the two results in Fig. 4, one may notice the merits and the 
weakness of the present approach for the quenched transition. The present approach ex- 
plains the periodic appearance of the large "shape factor" and well presents the period of 
its occurrence, but the details of the evolution is not exact. This discrepancy is related 
to the "patched" trial wavefunction (@]) and (|15|) at y = 3. We have chosen this artifi- 
cial patching method because of its simplicity. A better approach may be to include the 
excited states of (jH) without introducing the patching, (fT3j) . One of the excited states, 
{Q\ip2,t) = {y — l)~'^^'^ (g^y ~ (Q\'^'^)' is an orthonormal wavefunction to {Q\ip,t). Be- 
cause of the symmetry of the potential (P), the odd function of q cannot contributes to the 
evolution. One may try the variational method by using the following trial wavefunction: 



{Q\ij,t) = N-' [{Q\^,t)+z{Q\^2,t)] 



(23) 



where 2; is a complex valued function of time and the normalization factor is iV^ = 1 + 
This wavefunction naturally includes the regions with y > 3 due to the contribution of 
the excited state. In this sense, the appearance of large "shape factor" (> 3) is the signal 
for the contribution of excited states in the time evolution of the systems. Generally, the 
accuracy of the approximation (4) increases as the potential varies slowly. We applied the 
present method to the case of scalar 0^ field theory in Ref . P| , and it would be interesting 
to apply more realistic quantum mechanical systems in second order phase transition. 

Another point we need to speculate is the large N limit. It was shown that the large 
N wave-function satisfies jTU] 



1 92 



2iV2 dz"^ 



+ u{z,N) 



<l>(z,r). 



(24) 
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N 

where t = Nt, = J2 = and u{z, N) = ^^"jjifr^^ + K-^^ " 4?- The present 

k=l 

approximation for symmetric state is given by 



rfr(^, T\tdr - - u{z, AT) 1^, r). (25) 

This is the same as Eq. Q with the change of parameters t, Q, m, V{Q, t) — > r, z, N'^,u{z, N). 
With the use of the trial wave-function ^ in the expectation value of 2^n^ is given 
by 



21/ 



y 

y 



Y + l 



2gV _ 4yg^np . . 

m m m ' ^ ^ 



where the first three terms in Eq. ()26|) are 0(1), and in the large limit the quantum 
mechanical effects on the potential effectively vanishes as 0{1/N'^). In the absence 
of this quantum mechanical term, the equations of motion in the large limits for the 
present quartic exponential approximation with y = 1 is the same as that of the Gaussian 
approximation centered at ^ 7^ 0. Since the Gaussian approximation was proven to be 
the same as large approximation ^T], the present approximation is equivalent to the 
large A^ approximation for A^ — 00. 
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